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ABSTRACT 

The observed rotation curves of dark matter-dominated dwarf galaxies indicate low 
density cores, contrary to the predictions of CDM models. A possible solution of this 
problem involves stellar feedback. A strong baryonic wind driven by vigorous star 
formation can remove a large fraction of the gas, causing the dark matter to expand. 
Using both numerical and analytical techniques, we explore the maximum effect of the 
feedback with an instantaneous removal of the gaseous disk. The energy input depends 
on the compactness of the disk, hence the specific angular momentum of the disk. For 
the plausible cosmological parameters and a wide range of the disk angular momenta, 
the feedback is insufficient to destroy the central halo cusp, while the inner density 
is lowered only by a modest factor of 2 to 6. Any realistic modeling of the feedback 
would have even lesser impact on dark matter. We find that no star formation effect 
can resolve the problems of CDM cusps. 

Key words: galaxies: dwarf — galaxies: formation — galaxies: kinematics and dy- 
namics — dark matter 



1 INTRODUCTION 

Dwarf spheroidal galaxies are excellent systems for testing 
current theories of dark matter. The large mass-to-light ra- 
tios ind icate that ev en their centers are dominated by dark 
matter (Mateo 199^). The velocity dispersion profile of stars 



and the rotation curve of gas provide clean measures of the 
dynamical mass at all radii. Recent observations of the gas- 
rich dwarfs show a nearly solid-body rotatio n curve, which 
implies a finite core inside a few h undred pc ( van den Bosch 
et al. 200C| ; |de Blok et al. 200l| ). This is in confiict with 
the predictions of the cold dark matter (CDM) models. Nu- 
merical simulations invariably produce a d iverging, power- 



law density "cu sp" p oc r ^ with 7 = 1 ([Navarro, Frenk 



5OO(Mvir/lO^M0)^'^^ pc (Bullock et al. 2001a) 



fc White 19971 ) or 7 = 1. 5 ([Moore et al. 1999b| ) within 



The additional disagreements of the CDM models with 
observations include (1) the prediction of too many dwarf 
satellites within large halos ( [Klypin et al. 1999| ; Moore et 
al. 1999a); (2) the triaxiality of halos in conflict with the al- 
most spherical cores of som e clusters of galaxi es, as inferred 
from gravitational lensing ( Tyson et al. 1998| ); (3) the effi- 
cient transport of the baryon angula r momentum, leading 
to to o small disks of spiral galaxies (Navarro & Steinmetz 
200c). Of these, the central cusps seem to be the most severe 
discrepancy. The CDM models have been very successful in 
matching observations on scales larger than a Mpc, but the 
smaller-scale problems have lead to a recent search for al- 



terna tives to CDM. The variants i nclude warm dark matter 
(e.g., |Bodc, Ostrikcr fc Turok 200]|) and self-interacting dark 
matter dSpcrgcl fc Steinhardt 2*0001 ). 

We examine a possible astrophysical solution to the 
cusp problem, the effect of star formation feedback. The en- 
ergy released in supernova explosions may heat and ionize 
the surrounding gas. In the event of an extremely powerful 
burst of star formation, the heate d gas may leave the dwarf 



Dekel fc Silk 1986). In an 



galaxy in a form of fast wind (e.g.. 
idealization of this problem, a significant fraction of the gas 
may be removed from the dwarf halo on a timescale shorter 
than the dynamical time. We consider the reaction of the 
dark matter distribution to the sudden loss of baryonic mass 
in the center. Supermassive black holes can also have inter- 
esting effects, somewhat different from the gas disks, and 
these will be discussed elsewhere. 



2 COSMOLOGICAL CONTEXT 

We explore the maximum effect of stellar feedback on the 
dark matter profile, assuming (1) all cool gas within the dark 
halo can be removed; (2) the amount of gas locked in stars 
is negligible; (3) the removal of gas is instantaneous. 

The initial conditions for the dark matter halo and the 
baryonic disk are set using the cosmological co ncordance 
model w ith Q,m = 0.35, = 0.65, Vli = 0.02 ( [Bahcall et 
al. 1999), where h = 0.65 is the Hubble constant in units of 
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100 km Mpc~^. This gives the average baryon fraction 
fb = Q.b/^m ~ 0.13. Some of the baryons associated with 
the halo of mass Afvir will never cool and collapse towards 
the center, some may be heated and ionized early on and 
escape the shallow potential well before the starburst. But 
since we are interested in the maximum effect, we assume 
that all the associated baryons settle into a disk with the 
mass Md = fb M^ii. 

The size of the rotationally-supported disk, r^, is de- 
termined by the angular momentum of the baryons (e.g., 
Mo, Mao & White 1998). Cosmological simulations show 
that dark matter halos acquire a log-normal distribution of 
the dimensionless angular momentum, A = J\E\^^' ^ /GAl^J^ , 
with the median A = 0.05 and dispersion a\ — 0.5 (Barnes & 



Efstath ^ou 1987| [Efstathiou et al. 19881 ). If the gas had the 
same angular momentum as the halo, the disk size would 
scale with A {rd ~ A'^ r^ir if disk is fully self-gravitating, and 
rd ~ Arvir in the no self-gravity case). 

Initially, in isolated halos the gas is supposed to ac- 
quire the same angular momentum per unit mass as the 
dark matter. As a result of hierarchical merging, the distri- 
bution of dark matter can become highly non-homogeneous. 
Also, the cooling and fragmentation of gas into dense clouds 
can transfer orbital angular momentum to the dark matter 



via dyn amical friction and direct torques (Navarro & Benz 



199l |). \ n the end, the gaseous disk may have a significantly 
lower value of the parameter \b X. We mimic the effect 
of angular momentum redistribution by considering three 
model disks with different A;,. 

The dark matter distribution self-consistently adjusts 
to the formation of the disk in the middle of the h alo. This 
process can be described as adiaba t ic contraction ( Blumen 
thai et | al. 198^ ; ^lores et al. 1993|; ^Icanton et al. 1997| ) 



which we discuss at length in §3.1. Such condensation of 
both baryons and dark matter has a significant effect on 
the central density. Figure |^ demonstrates how adiabatic 
contraction increases the rotation curve in the inner galaxy 
over that predicted by cosmological simulations and inten- 
sifies the discrepancy with observations. 

After the gas blowout, the final distribution of dark 
matter will depend on the initial disk size in the following 
way. If the disk was initially relatively extended {rd ^Ts) the 
re-expansion phase is driven by the energy input, which is 
proportional to the binding energy of the disk, oc Md/rd- If 
the disk was very compact (r^ ^ rs) only the total removed 
mass Md affects the results. The new equilibrium profile af- 
ter the removal of the entire baryonic disk must be more 
extended than that prior to adiabatic contraction, because 
any rapid expansion increases the actions of dark matter 
orbits. The faster the expansion, the stronger is the effect. 
The energetics o f gaseous winds have been studied ex- 



tens ivel y in the past (Dekel & Silk 1986; 



van den Bosch et al 



200C ar d references therein). In general, dwarf galaxies have 



shallow potential wells and it is easy to create a strong wind 
after a few supernova explosions. We shall simply assume 
that the SNe ejecta energy is sufficient and is efficiently de- 
posited into the gas, to facilitate a rapid expulsion of nearly 
all the initial baryons with very little locked into stars and 
stellar remnants. 



> 




Figure 1. Rotation curve before and after the formation of bary- 
onic disk and adiabatic contraction of dark matter, for three val- 
ues of the angular momentum parameter, A = 0.05,0.01,0.0025. 

3 ANALYTICAL MODELS 

3.1 Two models of adiabatic contraction 

Wc consider the problem in a spherically-symmetric a pprox- 



imation. For the initial virialised halo model we use a Hern _ 
quist (1990| ) profile, which is a speci al case of a family of 



analytical double-power-law models ( Zhao 1996| ) with the 
mass Mdm{f'i) enclosed within radius ri given by 

Mdr^(ri) 



\rs +riJ 



3-7 



(l-/6)Afvir Ws+rJ ' ' ^' 

where is the scale radius, and 7 is the slope of the cusp. We 
adopt a concentration parameter rvir/^s = 10, appropriate 
for the dwarf halos at high redshift, but the results at small 
radii are insensitive to the actual value of this parameter. 
The pr ofile has the same inner slope ( 7 = 1) as the NFW 
model (Navarro, Frcnk & White 1997) and deviates from 



it only in the outer parts which are not interesting for the 
present problem. The Hernquist model is more convenient 
both for analytical modeling and for generating the initial 
condition of an isolated fixed-mass halo. 

During adiabatic contraction, circular orbits conserve 
the absolute value of the angular momentum, GM{r)r, and 
the enclosed dark matter mass, Mdmir). An orbit can be 
labeled interchangeably by its initial radius ri or by the 
specific angular momentum j, where j and ri are related 
by 



f ^r,GMdm{n){l- hY 



(2) 



Here we use the constant factor (1 — fb)~^ to include the 
mass of baryons assumed to have the same initial distribu- 
tion. The orbit with an angular momentum j has a post- 
contraction radius rj given by 



GMdm{r^) + GMbirj)' 



(3) 
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density of the disk would be too high and the scale length too 
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Figure 2. Baryon mass profile according to our angular momen- 
tum prescription (solid line, cf eq. |^ with 7 = 1 and n = 4/3), 
in comparison with the universal profile (dashed line) and the 
observed disk galaxies (symbols), where Jmedian is defined such 
that 50% of the disk mass has the specific angular mo rnentiim 
j < i„„ A:.,„. The observed disk pr ofiles are adopted from van den 
Bosch, jjurkcrt fc Swaters (2001 ) with the stellar mass-to-light 
ratio M/Lr = IMq/Lj^ q. 



where Mi,{r-j) is the mass of baryons enclosed within the 
contracted orbit, which we need to specify separately. 

We will use two independent models to evaluate adia- 
batic contraction. 

1) Using the observations of disk galaxies as a guide, it is 
common to assume that the disk would have an exponential 
surface density profile. The enclosed disk mass is 



Mi,,cxp(rj) = .fbA'/vir [1 - (1 + y) exp(-y)] , 



y 



where the scale length Vd is computed nu merically by match- 



ing the total angular momentum (e.g.. Mo, Mao & White 
199^ri^eir eq. 28). In this approach, r_, is an implicit func- 
tion of Ti or j, and the equation must be solved iteratively. 

2) Here we also take another approach, assuming the 
disk mass Adbixj) can be expressed as a function Mb{< j{ri)) 
of the initial specific angular momentum j{ri), where j{ri) 
is given by equation The advantage is that the r.h.s. of 
equation (^ can now be expressed as an explicit function of 
ri, and hence rj is computed directly without iteration as in 
the previous approach. 

Since the ^-component of the specific angular momen- 
tum originated from tidal torquing of neighbouring halos, 
hence it was distributed in the same way in baryons and 
dark matter, we could adopt the so-called universal angu- 
lar momentum profile, Mi,(< j) — /^Mvir j/(0.8j -|-0.2jmax), 
s een in dissipationless N-body simulations of ACDM models 
(Bullock et al. 2001h), where jmax is the maximum specific 
angular momentum determined by the spin parameter A. 
Unfortunately, if baryons with such a distribution collapse 
to a disk conserving angular momentum, then the central 



low compared to the observed disks ( Bullock et al. 2001b 

Swaters 200 1[ ). The angular mo 



van den Bosch, Burkert 



mentum distribution of such a disk is also very different from 
that of exponential disks, as shown in Fig. ^ This mismatch 
between the present day exponential disks and the simulated 
initial baryon angular momentum distribution is an intrigu- 
ing problem, which might argue for the transfer of angular 
momentum by bars or spiral waves in the disk or unvirialised 
lumps in the halo. We will not digress into this issue since 
the goal of our analytical and N-body modeling is to explore 
the maximum effect of blowing out the exponential disks. 

We have decided to make an alternative choice of the 
functional form of the initial cumulative distribution of the 
specific angular momentum of baryons. 



Mt{< j) 



= 1 



6-27 



(5) 



where 1.33 > n > 1.2 for a general cusp slope 1 < 7 < 1.5 
of the initial halo, and jo is a characteristic specific angu- 
lar momentum, related to the median angular momentum 
by jmcdian = Jo(ln2)^^". This parametrization can describe 
the angular momentum profile of observed disks fairly well 
(cf Fig. ^), albeit with significant scatter. It also creates an 
approximately exponential disk after adiabatic contraction. 
This can be understood roughly since j = rjVj oc rj in a 
model with a nearly fiat rotation curve, so that equation 
(^) corresponds to an almost exponential distribution in ra- 
dius Vj . Also, our angular momentum profile reproduces the 
asymptotic profile Mt oc j'*'^'^ cc r| near the center, as ex- 
pected for the exponential disk embedded in a Hernquist 
halo with 7 = 1 and n = 4/3. The characteristic specific 
angular momentum jo is determined by the integrated total 
angular momentum J of the initial halo. 



d[l-exp(-a;7)] J AGM^Y/l^l 



-1/2 



virJO 



JO 



,(6) 



where the l.h.s. can be reduced to the gamma function, r(l-|- 
1/n). Since the baryons have initially the same distribution 
as dark matter, the total energy of the system is 



\E\ = {i-hy 



dr 



dMdm GMd„ 



dr 2r 4(5 - 2'y)rs 
and we can express jo in terms of the A parameter, 
AVGMvirr. 



Jo = 



5 = r(l + l/n)/V4(5-27). 



(7) 



(8) 



At small radii the initial mass distribution of dark mat- 
ter is 

Mdrr.-^ (I- fb)M^ir{^y)' \ U ^ Q, (9) 

and the initial specific angular momentum 
j OC VMdnir^ oc r[ 



oc M 



(10) 

So with our choice of n, the asymptotic angular momentum 
distribution of the baryons scales the same way as that of 
the dark matter: 



A/d„ oc Ml 



b oc J 



0. 



(11) 



It follows then that the mass profiles of the baryons and 
dark matter are similar before and after the contraction. 



© RAS, MNRAS , 



4 Gnedin & Zhao 



3.2 Analytical models of slow wind and fast wind 



-p 0.6 




Figure 3. Baryon mass profile according to our angular momen- 
tum prescription, in comparison with the exponential disk profile 
(solid line; cf. eq. The radii are normalized to the scale length 
Tjj computed numerically for three cases: A = 0.05 (dots), A = 0.01 
(short dashes), A = 0.0025 (long dashes). 



Mh oc Mdm oc r. 
where 

2 

Tj OC j ^ — ' OC r. 



3 — 'Y -n 3 — "t* 

OC J OC ''j I 



(12) 



(13) 



is the post-contraction radius. The contraction factor is fi- 
nite at small radii, and the cusp slope of the darlc matter 
does not change. To summarize, at small radii after the con- 
traction we have 



Mb 



(1 - /6)Mvi, 



3-7 



Tj C ^Ti, 

where 



(14) 

(15) 
(16) 

(17) 

At a given radius rj , the enclosed dark mass and density are 
enhanced by a factor C'^ For an initial Hernquist halo, 
-y = 1, n — 4/3, and ^ = 0.265. For the disk parameters 
fb = 0.13, A = 0.05 — 0.01, the inner dark matter density is 
enhanced by a factor — 4 — 160. 

It is interesting to note that the model with 7=1 
produces an asymptotic mass distribution at small radii, 
Mt oc r|, similar to that of an exponential disk (cf eq. |^. 
Figure ra shows the computed mass profiles of baryons for a 
family of models with different Xb- Our analytical models re- 
produce exponential profiles at 10% level within 0.5 — 5 scale 
lengths. The difference in the disk potential is even less and 
therefore the adiabatic contraction of dark matter is very 
similar. For all three models, the dark matter profiles com- 
puted with the two methods are essentially indistinguishable 
between 10~^ and 10 rg. 



The expulsion of gas can be worked out analytically in two 
Umiting cases: the adiabatic and instantaneous regimes. We 
shall refer to these as the slow wind and the fast wind, re- 
spectively. The outfiow can play a significant role in massive 
galaxies as well as in dwarf galaxies. Our analytical descrip- 
tion is similar to that used in estimating the maximum wind- 
induc ed expansion of the dark halo of the Milky Way (Zhao 



2002) 



To gain more analytical insight, we first consider the 
slow wind. In this case, a circular orbit of radius Vj expands 
to a circular orbit of radius r/ while conserving the spe- 
cific angular momentum -y/ GM{r)r. So the expansion factor 
is inversely proportional to the combined enclosed mass of 
baryons and dark matter, i.e., 



Mam{rj) + Mb{rj) - M^{rj) 



(18) 



where Mbirj) and Mdm{rj) are the baryon and dark mat- 
ter masses within radius rj, and Mw(rj) is the amount of 

Alternatively, we can write 

(19) 
(20) 



baryons "gone with the wind" 
where 



Md^{rj)+Mb{rj)' 

measures the blown-away mass as a fraction of the dynam- 
ical mass. The expansion factor rf/r^ is larger in the inner 
region where the baryons are concentrated and is minimal 
at the virial radius, where we have only a modest change, 



= Sft < 0.1, 



(21) 



where < (5 < 1 is the total fraction of baryons blown away. 

Although we do not know of any close form for the effect 
of fast wind, we could modify equation (|l^) for slow wind to 
mimic the effects of wind of arbitrary strength. We suggest 
the following simple formula 



< it < 1, 



(22) 



[l-.Urj)]-il-2'=Sfb)- 

where fc is a tunable parameter to model the rapidity of the 
wind. For example, slow wind is a special case with A: = 
(cf eq. p^ . The instantaneous wind is modeled with k — 1, 
following the rule of thumb that a system becomes unbound 
if it loses half of the mass instantaneously. To verify this, we 
substitute k = 1 and fm{rvir) — > 1/2 in equation (p^ ) and 
find that the r.h.s. factor 1 - 2''6fb = 1 - 2/„(rvir) ^ 0, 
hence r/ — > oo. For the winds of moderate speed with A; < 1 
and moderate mass-loss with /u,(rvir) < 2"*^, the system 
remains bound and expands by a certain factor at large 
radii. At small radii, where the baryons always dominate 
with fuj{rj) — > 1, the particles may not necessarily escape 
to infinity because they slow down when they cross the or- 
bits of outer particles. As we will show in §^ this empirical 
parametrization turns out to approximate well the results of 
numerical simulations. 

To work out the properties of the model at very small 
radii, we substitute equations ( [l4| - [l^ ) in equation ( pc| ) and 
obtain 
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U^C^~"SC-'h. (23) 

Substituting equations ( [l^ , p^ in equation (p^, we find that 
in the center the final post-wind radius r/ is related to the 
initial pre-contraction radius by 



1 + (^A- 



[1 + (eA-"(l -S)- 1)M'+- (1 - 2Hft)- 



.(24) 



The slope of the cusp remains the same. By comparing the 
dark matter mass within a given radius, we find that the 
density is reduced by a factor F'^~'^ . 



3.3 Maximum effect on the final halo density 

The strongest effect on the halo is achieved in the limit of 
complete removal of the baryonic disk, i.e., 5 ^ 1. In this 
case we have for all radii 



(1 - .h) 



1 + Mb/Ma„ 
l-2Vi> 



(25) 



Interestingly, if the disk forms slowly and is removed slowly 
as well (i.e., with k — 0), the halo expands slightly by a 
constant factor 



— = {1 — ft) ^ — const, 



(26) 



as expected from adiabatic invariance. The result is the same 
if all the baryons were removed from the initial halo distri- 
bution without going through the phase of disk formation. 

The effect is maximized when a very dense disk (A ^ 
0.05) is removed instantaneously [k = 1). If the initial halo 
is a Hernquist model (i.e., 7 = 1, n = |, and ^ = 0.265), 
then close to the center the halo will expand by a factor (cf 
eq.pi) 



F ' 



1 



(l-/6)(l-2/5) 



0.265 \ 3 
A ) 



(27) 



relative to the original halo model. Here we assume a uni- 
versal baryon fraction ft = 0.13 to make a massive disk. To 
summarize, the density of the innermost dark matter can 
drop at most by a factor = 2 — 6 for A = 0.05 — 0.01. 



4 SPHERICAL SHELL SIMULATIONS 

To confirm the analytical results, we compute the dynamics 
of dark matter using the spherical shell method of Gnedin & 



□striker (1999). The algorithm, originally due to L. Spitzer 



and collaborators, utilizes spherical geometry to achieve 
high resolution in the center. Dark matter particles are dis- 
tributed on concentric spherical shells at radii with the 
mean energy Ek and angular momentum jk- Shells move 
in the radial direction with the velocity Vr,k and can freely 
cross each other. 

We have extended the code of 



Gnedin & Ostriker (199£) 



by including variable shell masses, mk, and improving the 
corrections for shell crossing. The code solves the equations 
of motion using the first integral, Ck'- 



„,2 I jk 



2GM{rk) 



(28) 




Figure 4. Dark matter density profiles for Model A (A = 0.05, 
/i, = 0.13, 5 = 1) before and after the gas removal. Dashed line is 
the initial Hernquist model, solid line is the halo profile after the 
adiabatic contraction. Symbols show a new equilibrium profile af- 
ter the removal of the baryonic disk, computed with 10^ spherical 
shells. Errorbars are the sampling uncertainty. Dotted line is the 
analytical prediction, eq. (^). 



where M{rk) is the mass enclosed within shell k, including 
a half of its own mass (accounting for the shell self-gravity): 



fc-i 



(29) 



=1 



if the shells are ordered in radius. The integrals Ck determine 
the turn-around points where the direction of the shell's 
velocity is reversed. At each time step, the code finds the 
locations of all shell crossings and adjusts the integrals Ck 
according to the changed enclosed mass. Because of such cor- 
rections, two-body interactions are effectively reduced and 
there is no need for force softening. As a result, the central 
region of the system can be probed with very high accuracy. 

In order to investigate the effect of angular momentum 
loss, we have run three models: Model A with the mean value 
of A = 0.05, and Models B and C where A deviates from its 
mean by 3cta (A = 0.01, Model B) and 6cta (A = 0.0025, 
Model C). 

We use the units G = 1, Mvir = 1, = 1, and there- 
fore the results can be scaled to any halo parameters. Each 
model is run with iV = 10*, 4 X 10'', and lO'^ shells to check 
numerical convergence. After the disk is removed, the mod- 
els are run for 40 half-mass dynamical times (tdyn = 8 in 
code units) to ensure that a new dynamical equilibrium is 
achieved. 

Figure ^ shows the initial halo profile and the new 
equilibrium distribution after the instantaneous removal of 
the whole disk in Model A. The scale length of the disk, 
rd = 0.085 rs, is already a small fraction of the scale radius 
of the halo. As expected, there is no change in the profile 
outside the cusp, at r > rs. But even inside the scale radius. 
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Q, 




10-2 IQ-i I 10 



r/r 



ure |6|) . The inner profile can be fitted witli an approximate 
isotliermal spfiere, p oc (r^ + r^)~^, wfiere Tc ~ 0.2 r^. Tliis 
core radius is significantly larger than the tiny initial scale- 
length of the disk, = 0.00041 r-^. However, it seems that 
the energy input has already saturated in Model C (the 
change in the binding energy is similar to Model B) and 
further contraction of the disk would not lead to a larger 
core. 

After the disk is removed, the dark matter particles with 
a large enough kinetic energy become unbound. At the end 
of the simulation, the halo in Model A loses 6.5% of its mass, 
while halos B and C lose 14% and 16%, respectively. (About 
half of those particles become unbound during the subse- 
quent evolution of the halo potential towards a new equi- 
librium.) However, as the disk is made progressively more 
compact, the heating caused by the disk removal is shared 
by a small number of particles in the center. They reach 
very high velocities and leave the system, whereas most of 
the particles remain unaffected. In the limit of an infinitely 
small disk, the final equilibrium profile would still be close 
to our extreme Model C. 



Figure 5. Density profiles for Model B (A = 0.01). Lines are as 
in Figure ^ 
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Figure 6. Density profiles for Model C (A = 0.0025). Lines are 
as in Figure ^. 



the new profile is fairly close to the initial model prior to adi- 
abatic contraction. The central density is affected at most 
by a factor of 2, in agreement with the analytical prediction 
(eq.H). 

Figure shows the density profile for Model B. The ef- 
fect is stronger but there is still no indication of a core form- 
ing instead of the cusp. The asymptotic solution increases 
towards the center as and describes the numerical result 
reasonably well. 

Finally, Model C may have achieved a central core (Fig- 



4.1 Artificial numerical effects 

As any N-body simulation, our models are affected to some 
extent by the artificial numerical relaxation. This relaxation 
arises from the small-scale particle encounters an d from the 
Poisson potentia l fiuctuations due to discreteness (Hernquist 
fc Barnes l"99o|; [Weinberg 199^). Even though the shell- 



crossing corrections in our code remove the small-scale er- 
rors, the large-scale fiuctuation modes are unavoidable. 

The relaxation time over which the particle energy 
changes by the order itself can b e estimated at the hal f-mass 
radius of the system, Rh, from dSpitzer fc Hart 197l[ ): 

N /"3XV2 



0.14 



In A \GM 



(30) 



where In A ~ InA'^ is the usual Coulomb logarithm. In the 
code units, G = 1, M = 1, and with Rh ~ rs = 1, we have 

U^O.l^. (31) 

This expression can be compared with the relaxation 
rate ( inverse of the relaxation time) derived by Weinberg 
(1993). His Fig. 3 expresses tr as a function of the system 
size divided by the Jeans length. For a virialised Hernquist 
model, this ratio can be 0.6—0.8 and reading of the plot gives 
tr ~ 25 — 50. It agrees wi th equation (|3l|) for the particle 
number A'' = 4000 used by [Weinberg (1993| ) . 

Thus, the large-scale potential fiuctuations can only be 
reduced by increasing the number of particles in the sim- 
ulation. We estimate that our results with A'' — 10^ shells 
are unaffected by the relaxation at the scale radius for 
t <^ tr ~ 800 « lOOtdyn. The rate of relaxation does not 
increase much in the inner parts, because it is dominated by 
the global modes and not by local encounters. 

Another numerical effect, which is more important for 
the inner shells than the outer shells, is the kinematic error 
due to a finite time step. We chose At = tdyn/N, so that each 
shell crosses on average one other shell. Since the density at 
r < goes approximately as p(r) oc r~^, the orbital period 
of a shell torhir) ~ tdyn (r/rs)^^^. This shows that if 10 time 
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steps are necessary to properly resolve an orbit, our results 



are valid down to a tiny radius r ~ (10/A'^) 



4.2 Comparison with previous work 



10" 



iNavarro, Hke & t'renJt (iyyti) studied the effect of instanta- 
neous disk removal using a three-dimensional N-body code. 
They also used a Hernquist model for the initial dark matter 
profile and imposed an external potential of the exponential 
disk. Instead of analytically calculating the adiabatic con- 
traction, they let the dark matter halo to adjust dynami- 
cally to the disk potential. Then, they removed the disk and 
studied the subsequent expansion of the halo. Navarro et al. 
found that the core did develop in the inner regions of all of 
their models. They fit a non-singular isothermal sphere to 
the inner 25% of the mass and found the core radii satisfying 
the following relation: Vc = 0.11 {Md/rd)^^^ rs, where Md is 
the mass of the disk in units of the halo mass (or fb/{l — ft), 
in our notation). 

The expansion effect that Navarro et al. find is much 
stronger than what our simulations or analytical modeling 
show. For example, their fit for the core radius would give for 
our Model C, Vc « 2.1, about ten times the value we find. In 
Models A and B the core radii should also have been easily 
detectable, according to the fit. 

We identify two possible reasons for the disagreement, 
force resolution and numerical rela xation effects. Firstly, a 
force softening at 0.03 rs, used by Navarro, Eke & Frenk 



(1996 ), could prevent an accurate calculation of the dy- 



namics of inner particles, especially in cases of very small 
disks. This would naturally cause a development of a core 
in the particle distribution, although the core radii found by 
Navarro et al. are larger than the softening length. Secondly, 
numerical effects might be important, as their simulations 



use only A'^ = 10 particles. The relaxation time, tr 



100 



time units, is 4 to 6 times shorter than the duration of their 
experiments. Even though Navarro et al. demonstrate that 
the numerical profile before the disk removal is stable, it 
may be more susceptible to the relaxation afterwards, when 
a smaller number of particles is left in the center. This may 
also lead to a simple statistical undersampling of the density 
distribution, posing as a 'core'. 

It is likely that the difference in the results is not re- 
lated to the geometry, ID vs 3D simulations. In our one- 
dimensional models, all shells conserve their initial angular 
momentum and do not experience any non-radial perturba- 
tions. But since the dark matter halo is not centrifugally sup- 
ported, precise conservation of angular momentum should 
not seriously affect the dynamics. Also, the remaining stel- 
lar distribution in dwarf spheroidal galaxies is observed to 
be roughly spherical, so we can expect the dark matter halo 
to be spherical as well, lending support to our ID models. 
Finally, we have checked that our different implementations 
of the adiabatic contraction gave almost identical initial con- 
ditions and therefore cannot be responsible for the disagree- 
ment. 

Note, that our simulations achieve much higher reso- 
lution, with ten times more particles and no force soften- 
ing in the center. But the statistical error-bars in Figure ^ 
inversely proportional to the square root of the number of 
particles in the bin, show that even our models cannot probe 
regions smaller than 0.01 — 0.03 rs. At larger radii, except 
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Figure 7. R otation curve of DDO 15 4 from a high-resolution 
HI data (from Carignan & Purton 199^ ) compared with the final 



circular velocities of Models A (dots), B (short dashes), and C 
(long dashes). The simulated rotation curves are normalized to 
pass through the maximum of the observed one, with no formal 
fitting attempted. 



for the extreme model C, we do not find any significant flat- 
tening of the dark matter distribution. 

Any realistic mod eling of stellar feedback wou ld have 



an even weaker effect. Lia, Carraro & Salucci (2000), using 



a more complex tree-SPH simulation of star formation in a 
dwarf galaxy, find that the central cusp remains intact. 



4.3 Application to DDO 154 

Our results demonstrate that stellar feedback is insufficient 
to reduce the central dark matter density. We illustrate this 
point for a prototypical low surface brightness galaxy DDO 
154. Figure Q shows the observed rotation curve of DDO 
154, which peaks at 48 km s~^. As this is perhaps the most 
robustly determined quantity, we normalize the simulation 
results to have the same value and location of the circular 
velocity. The simulated curves correspond to the final equi- 
librium density profiles of Models A, B, and C. They are 
clearly more concentrated than allowed by the data. 

In addition, the scale radii of the initial halo in Figure ^ 
increase as the effect of feedback gets weaker (corresponding 
to the larger disk). For Model C rs — 2.9 kpc, but for the 
less perturbed models B and A it is ra = 3.6 and 5.6 kpc, 
respectively. The virial halo mass is about the same, around 
1.3 X 10^° Mq. This m ass is lower than the best-fitt ing value 



(~ 3 X 10^" Mq) from [van den Bosch et al. (2000[ ), and yet 



the resulting rotation curve is well above observations. 

Note, that we do not attempt any accurate fitting of the 
rotation curve. The beam-smearing effects are likely to be 
unimportan t (beam size is less than 1 kpc; van den Bosch 
et al. 200C), and therefore the small error-bars in the in- 



ner parts place very strong constraints on the allowed dark 
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matter profiles. Selato fc Sommer-Larsen (199S: ) used a set 
of N-body simulations to address this issue. They selected 
three halos with M ~ 10^" Mq at 2: = from a random 
realization of a 4/i~^ Mpc box with the SCDM power spec- 
trum. Varying the ar nount of gas and the concentra tion of 
the dark matter halo, [Gelato fc Sommer-Larsen (1999 ) could 
find a satisfactory fit to the inner rotation curve of DDO 
154 only (1) if the present disk contained 3 times more gas 
than is observed and (2) if the halo had a low concentration, 
rvii/rs — 4. However, none of their models could fit simul- 
taneously the inner and outer parts of the rotation curve. 



5 SUMMARY 

We have explored the effect of maximum feedback on the 
central density of dark matter halos of the gas-rich dwarfs, 
combining analytical limits with numerical simulations. The 
expansion of dark matter after the removal of the entire disk 
is controlled by the mass and compactness of the disk. For 
a wide range of the baryon angular momenta, we find the 
effect to be modest, at most a factor of 2 to 6 reduction in 
the central halo density. This is hardly enough to bring the 
models into agreement with the observed solid-body rotation 
curves, as we demonstrate for the case of DDO 154. 

We conclude that the slowly rising rotation curves are 
likely to be a genuine problem of CDM models. It might 
be necessary to consider other possible solutions, such as 
the effect of merging of massive black holes or an unusual 
property of dark matter particles. 
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